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ABSTRACT 

The recently discovered magnetic Herbig Ae and Be stars may provide qualitatively 
new information about the formation and evolution of magnetic Ap and Bp stars. We 
have performed a detailed investigation of one particularly interesting binary system 
with a Herbig Ae secondary and a late B-type primary possessing a strong, globally 
ordered magnetic field. Twenty high-resolution Stokes V spectra of the system were 
obtained with the ESPaDOnS instrument mounted on the CFHT. In these observa- 
tions we see clear evidence for a magnetic field in the primary, but no evidence for 
a magnetic held in the secondary. A detailed abundance analysis was performed for 
both stars, revealing strong chemical peculiarities in the primary and normal chem- 
ical abundances in the secondary. The primary is strongly overabundant in Si, Cr, 
and other iron-peak elements, as well as Nd, and underabundant in He. The primary 
therefore appears to be a very young Bp star. In this context, line profile variations 
of the primary suggest non-uniform lateral distributions of surface abundances. In- 
terpreting the 0.63995 ± 0.00009 day variation period of the Stokes / and V profiles 
as the rotational period of the star, we have modeled the magnetic held geometry 
and the surface abundance distributions of Si, Ti, Cr and Fe using Magnetic Doppler 
Imaging. We derive a dipolar geometry of the surface magnetic field, with a polar 
strength B d = 1230 G and an obliquity = 57°. The distributions Ti, Cr and Fe 
are all qualitatively similar, with an elongated patch of enhanced abundance situated 
near the positive magnetic pole. The Si distribution is somewhat different, and its 
relationship to the magnetic field geometry less clear. 

Key words: stars: magnetic fields, stars: abundances, stars: chemically peculiar, 
stars: evolution, stars: individual: HD 72106 



1 INTRODUCTION 



Strong, globally organised magnetic fields have recently 
been reported in a few Herbig Ae and Be (HAeBe) stars 



dDonati et al.l Il997 


: iHubris et al.l 12004; 


Wade et all 120051, 


120071: Catala et al. 


120071; Alecian et all 


2008a bh. HAeBe 



stars are pre-main sequence stars of intermediate-mass 



which evolve to become main sequence A and B stars. 
HAeBe stars have A or B spectral classes, display emission 
lines and infra red excesses, and a re usually found with nebu- 
losity nearby l|Vieira et alj|2003h . The detection of magnetic 
fields in HAeBe stars is of particular interest as it hints at 
a connection to the main sequence magnetic, chemically pe- 
culiar Ap and Bp stars. 
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Ap and Bp stars display strong magnetic fields, with 
typical strengths of a few hundred gauss up to a few tens of 
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kilogauss, and globally ordered geometries of approximately 
dipolar topology. These stars also display strong and distinc- 
tive chemical peculiarities, particularly overabundances of Si 
and iron peak elements, occasionally in excess of 2 dex, and 
even greater overabundances of some rare earth elements. 

The source of the magnetic fields observed in Ap/Bp 
stars is not well understood. The two major competing field 
origin theories propose, on the one hand, that the field is a 
relic of an earlier stage of stellar evolution, now frozen into 
the plasma of the star or, on the other hand, that the field is 
generated contemporaneously by a dynamo. Also unknown 
are the details of the formation of the observed chemical 
peculiarities. Chemical peculiarities are believed to be the 
result of atomic diffusion leading to a chemically st ratified 
stellar envelope l|Michaudlll970l ; iMichaud et al.lll98ll ). How- 
ever, the details of this process, particularly for individual 
elements, and t he impact of magneti c fields are not fully 
understood (e.g. lAlecian fc StiftJ 120071 ). In this context it is 
of great interest to identify and characterize the evolution- 
ary progenitors of Ap/Bp stars, as these objects can provide 
critical information on both the origin of magnetic fields and 
the formation of chemical peculiarities. 

The recently discovered magnetic HAeBe stars have 
been propose d to be pre-main sequence progenitors of 
Ap/Bp stars ( Wade et al.l l2005h . Thus a detailed investi- 
gation of these stars can shed light on these questions, as 
well as potentially provide some information about magnetic 
braking and magnetospheric accretion in intermediate-mass 
stars. 



Canada France Hawaii Telescope. A longitudinal field of 
195 ± 40 G (i.e. 4.9cr) was deduced from the FORS1 spec- 
trum. No significant magnetic field was detected in the sec- 
ondary (65 ± 55 G was observed). Although the average 
longitudinal field in the ESPaDOnS spectrum of the pri- 
mary was consistent with zero, a clear circular polarization 
signature was detected with a very high degree of signifi- 
cance, unambiguously indicating the presence of a surface 
magnetic field. 

A more conservative re-analysis of the same FORS1 
data by IWade et al.l (|2007l ) did not confirm t he longitudi- 
nal m agnetic field detection at 3cr confidence. Wad e et al.l 
(120071 1 reported longitudinal fields in HD 72106A of 166 ±70 
G from Balmer lines and —11 ± 91 G from metallic lines, 
and fields in HD 72106B of 52 ±90 G from Balmer lines and 
3 ± 122 G from metallic lines. However they stressed that 
ESPaDOnS observations show a clear circular polarization 
signal, implying that HD 72106A is definitely magnetic. 

Due to its young age, probable binarity, and magnetic 
properties, the HD 72106 system is a compelling target for 
further study. While other binary systems containing a mag- 
netic intermediate-mass star are known, the apparently very 
young age and wide visual separation of HD 72106 make this 
system uniquely interesting. We therefore obtained 20 high- 
resolution spectropolarimetric (Stokes / and V) observations 
of the system, with the aim of studying the system's binarity, 
rotational properties, magnetic field, chemical abundances, 
and surface abundance distributions. 



2 HD 72106 



HD 72106 is a double star system in the constellation Vela, 
with a 0.805 arcsec separation between components (|ESA| 
Il997j ). The brighter star, HD 72 106A ('the primary '), was 



le p rii 

identified as a magnetic star by I Wade et al.l (|2005l ) . while 
the fainter star, HD 72106 B ('the secondary') , was iden- 
tified as a HAeBe star by IVieira et all |2003|). Th e HD 
72106 system was observed by Hipparcos |ESA| [i997). and 
included in the rece nt re-reduction of Hipparcos data by 
Ivan Leeuwenl l|2007al fbT) who found a parallax of 3.60 ± 1.14 
mas, placing the system at a heliocentric distance of 278lgy 9 
pc. The system was identified as having an infrared ex- 
cess, bas ed on IRAS (the I nfrare d Astronomica l Sate llite) 
data, by lOudmaiier et alj (|l992h . iTorres et all (|l995l ) ob- 
served HD 72106 and noted the presence of emission in the 
Ha Balmer l i ne of the combined spectrum of the system. 
ISchiitz et all (2005) obtained infrared spectroscopy of HD 
72106 system and examined circumstellar abundances of a 
number of dust g rain s pecies. 



IVieira et alj (|2003h associated the system with the Gum 
Nebula star forming region. They observed weak Ha emis- 
sion and a small contribution from dust in the spectral en- 
ergy distribution of HD 72106B. This, they hypothesized, 
was due to HD 72106B being an evolved HAeBe star that 
has cleared most of its circumstellar envelope. 

The dis c overy of a magnetic field in in HD 72106A by 
IWade et all {2005) was based on spectropolarimetry from 
FORS1 (FOcal Reducer/low dispersion Spectrograph) at the 
Very Large Telescope and from ESPaDOnS (Echelle Spec- 
troPolarimetric Device for the Observation of Stars) at the 



3 OBSERVATIONS 

Observations used in this paper were obtained with the 
ESPaDOnS instrument on the 3.6 m Canada- France-Hawaii 
Telescope (CFHT). ESPaDOnSSis a high resolution echelle 
spectropolarimeter, with a resolving power of about 65000 
in spectropolarimetric mode, and nearly continuous wave- 
length coverage from 370 to 1050 nm. This instrument uses 
a fiber fed design, with a polarimeter module mounted at the 
Cassegrain focus of the CFHT and an echelle spectrograph 
located in the Coude room. 

Twenty Stokes / and V spectra were obtained over a 
period of two years, and are summarized in Table Q] "Nor- 
mal" readout in spectropolarimetric mode was used, and the 
atmospheric dispersion corrector was employed for all ob- 
servations. Magnetic and non-magnetic standard stars were 
observed as part of this campaign, and produced results in 
excellent agreement with the literature. 

With careful guiding and monitoring, observations of 
the individual components of HD 72106 were acquired on 
a few nights with particularly good seeing. This was chal- 
lenging, however, as the double star system has a separation 
of only 0.805 arcsec and ESPaDOnS has a pinhole with a 
diameter of 1.6 arcsec. On most nights the atmospheric con- 
ditions were not good enough for us to resolve the individ- 
ual components of the double star system. In these cases the 
photocenter of the system was observed, providing combined 
spectra of the system. 



^ See http:/ /www. cfht.hawaii.edu /Instrumcnts/Spcctroscopy/Espadons/ 

for details. 
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Table 1. Log of observations of the HD 72106 syst em obtained with E SPaDOnS. Integration times refer to the complete set of 4 sub- 
exposures used to produce each Stokes V spectrum l|Donati et al]|l997h . Peak signal-to-noise ratios are reported for the original reduced 
spectra with 1.8 kms" 1 spectral pixels. Magnetic measurements for combined observations refer to results for the reconstructed spectra 
of the primary. The letter D in the LSD detection column indicates that the Stokes V LSD profile is inconsistent with the null field 
hypothesis at the 99.999% confidence limit, indicating the presence of a magnetic field. The letter N indicates that the inconsistency was 
at the 99.9% limit or less, and hence cannot be take as strong evidence for the detection of a magnetic field. The letters MD indicate 
a 'marginal' detection, w ith a value between 99.9% and 99.999%, and are considered sugge stive of a ma gnetic field, but not conclusive. 
l|Donati et al.lll99l \l99% . Note that the observation from Feb. 2005 was first published by I Wade et~aH (2005h . 
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The observatio ns were reduced with Libre-ESpRIT 
l|Donati et al.lll997l . and in preparation). This is a nearly- 
automatic dedicated data reduction package for ESPaDOnS, 
which performs a complete calibration and optimal spec- 
trum extraction. Continuum normalization of the resulting 
unnormalized spectra was performed manually, first by com- 
puting a running average on each order of an observed spec- 
trum, then selecting points in the continuum (i.e. maxima) 
of the smoothed order, and fitting a polynomial (typically of 
degree 4 or 5) through the selected points. The original ob- 
servation was then divided by the polynomial to produce a 
normalized spectral order, and the process was repeated for 
each order in the echelle spectrum. Similarly, Stokes V spec- 
tra were normalized by the Stokes I continuum polynomial, 
producing a V/I c spectrum. 



4 SPECTRUM RECONSTRUCTION 
4.1 Procedure 

The two individual spectra of HD 72106A show strong 
metallic lines with significant variability. In the two obser- 
vations of HD 72106B on its own we see no metallic absorp- 
tion line variability, down to the level of the noise in our 
observations. However, we do see emission in two lines of 
HD 72106B: clear, variable, emission is present in Ha, and a 
small amount of emission in the O I 7771 A triplet, as shown 
in Figure [1] Thus we confirm the spectroscopic properties of 
the secondary leading to its classification as a HAeBe star. 



The emission in HD 72106B varies detectably from night 
to night. Stronger variations are observed on timescales of 
months to years. No detectable change is observed between 
observations obtained on the same night. A period analysis 
of the emission line variations yields no significant periodic- 
ity in Ha. 

As a consequence of the ~0.8 arcsec angular separation 
between the components of HD 72106, the majority of our 
spectra of HD 72106A & B are of the combined light from the 
system rather than light from the individual components. 
The analysis of observed spectra is much more tractable 
when one is dealing with light from only one object. Thus 
we attempted to reconstruct the spectrum of one star from 
the combined spectra of the system. 

As the spectrum of the secondary does not appear to 
vary (outside of emission lines), we can subtract it from 
the combined spectra and reconstruct the spectrum of the 
variable primary. We begin by modeling the normalized flux 
I\ t T a t any point in the combined spectrum according to: 



I\,tLt — Ix.iLi + IX.2L2 



(1) 



where 7a, 1 and 7a, 2 refer to the normalized flux spectrum 
of the primary and secondary respectively at wavelength A, 
and Li and L2 refer to the luminosity of the primary and 
secondary respectively. Thus the product Ix,iLi gives the 
total flux of the primary. 7a, t and Lt are the observed nor- 
malized flux spectrum of the system and the total luminosity 
of the system respectively. The total luminosity can be writ- 
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Figure 1. Emission and variability in the Ha Balmer line and the O I 7771 A triplet of HD 72106B (the HAeBe secondary), shown for 
the two single spectra of HD 72106B. The smooth line in the O I 7771 A panel is an LTE synthetic profile calculated assuming solar 
abundances. Strong, variable emission can be seen in Ha, while mostly we see only emission infilling in O I 7771 A. No other lines in the 
spectrum of the secondary show any evidence of variability. 



ten as Lt = L\ + £2- It is elementary to rearrange Eq. |T} 
to obtain the normalized spectrum of the primary, Ix.i- 



(2) 



Thus the normalized spectrum of the primary can be re- 
constructed from an observed spectrum of the system (Ix,t) 
and a spectrum of the secondary on its own (7>, 2)- Formally 
the luminosity ratio should also depend on wavelength, how- 
ever since the two compon ents differ in temp erature by only 
2000 K (see SectionE21 or lWade et"aLll2007l ) the wavelength 
dependence is not strong, and thus has been neglected. Nev- 
ertheless, the magnitudes used to calculate the luminos- 
ity ratio must approximately correspond to the wavelength 
range of interest in the observed spectrum. The magnitude 
difference between components (0.62 ± 0.02 in the Tyc ho V 
band) is well known from Hippa rcos and Tycho data (|ESA| 
ll997l ; lFabricius fc Makarovll 20001) . thus the ratio of luminosi- 
ties at ~ 5000 A is well established. 



4.2 Quality Control 

To investigate the importance of the assumption of a 
wavelength-independent luminosity ratio, we used ATLAS9 
flux distributions corresponding to the atmospheric param- 
eters of the primary and secondary stars (Table [5}. The 
synthetic flux distributions show that the monochromatic 
luminosity ratio varies by 15% over the ESPaDOnS spectral 
range. However, since the reconstructed spectrum is essen- 
tially a difference spectrum, variations in the luminosity ra- 
tio correspond to variations in line depth Thus a typical line 
with a depth of 10% of the continuum, at the red or blue 
end of the spectrum, would have an inaccuracy of ~1.5% of 
the continuum, which is near the noise level in the obser- 
vations. For most of the lines used in this study this sys- 
tematic error is even smaller. Comparisons were performed 
using subtracted spectra calculated assuming different lu- 
minosity ratios, as well as comparisons of our reconstructed 
spectra to observations of the individual components. These 
showed that systematic errors due to our assumption of a 



wavelength independent luminosity ratio generally are well 
below the noise level in our spectra (with the exception of 
the Balmer lines, which were not included in our analyses). 

Stokes V spectra of the primary were reconstructed by 
attributing the observed polarization signature to the pri- 
mary and scaling the spectrum based the luminosity ratio 
of the stars. The assumption that the polarization signal 
is due solely to the primary is quite reasonable as we see 
no sign of circular polarization in our spectra of just the 
secondary. For the purposes of the subtraction, the higher 
S/N observation of the secondary (HJD 2454164.84650) was 
used. 

Contamination of the individual spectra of the sec- 
ondary by stray light from the primary is potentially a se- 
rious source of error in our analysis. However, in our obser- 
vations this contamination appears to be minimal. Obser- 
vations of the stars obtained on two different nights, under 
different seeing conditions, produce identical line profiles (in 
stable regions of the spectrum) at the level of the noise. If 
there were a significant amount of contamination, one would 
expect the contamination to vary with the seeing conditions. 
Thus the fact that we do not see any such variation suggests 
that the spectra are not detectably contaminated. A theoret- 
ical estimate of the fractional contamination in our spectra, 
based on the seeing and geometry of the observations, pro- 
vides further support to this claim. The stars were observed 
with the center of the target star on the edge of the pin- 
hole (0.8 arcsec from the pinhole center) and the other star 
placed outside the pinhole, ~1.6 arcsec away from the pin- 
hole center. Assuming a Gaussian point spread function for 
both stars, with a full width at half maximum equal to the 
seeing, one can integrate through the pinhole to estimate 
the flux from each star in our observation. Normalizing the 
Gaussians by the relative luminosity of the stars, and using 
the worst seeing conditions for our observations (0.9 arcsec 
for the primary and 0.7 arcsec for the secondary), we ob- 
tain a theoretical worst-case contamination of 6% (of the 
total observed flux) for the observation of the primary and 
7% (of the total observed flux) for the observation of the 
secondary. For contamination by a typical strong line (10% 
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Figure 2. Comparison of an observation of just HD 72106A 
(black, 5 March 2007), a combined observation of HD 72106A & 
B (dark gray/blue, 11 January 2006), and the corresponding re- 
constructed spectrum of HD 72106A (light gray/green) obtained 
using the procedure described in Section [4] A good match of the 
reconstructed spectrum to the observation of the primary can be 
seen in the core (upper left) and wings (upper right) of H/3, as 
well as in weakly variable metallic lines (lower panel). 

of the continuum), a deviation of roughly 0.6% of the con- 
tinuum would be seen in the observation, which at about 
the noise level of our spectra. Thus we conclude that such 
contamination is not a serious problem. 

Calculation of the spectrum of the primary from Eq. (J2J) 
was performed for all our combined observations of the sys- 
tem, allowing us to reconstruct a sequence of 18 spectra 
of the primary, typical with a peak S/N ratio of 140 (at 
5200 A). Comparisons of the reconstructed spectra to ob- 
servations of the primary on its own show excellent concor- 
dance in the Balmer lines, and in weakly variable metallic 
lines, as shown in Figure [5] Furthermore, this procedure can 
be applied to recover the spectrum of the secondary (with 
a typical peak S/N of 90) if we use regions containing only 
weakly variable lines of the primary. Such experiments also 
yield good agreement with observations of the secondary on 
its own. Thus we conclude that the spectrum reconstruction 
is robust, and yields high-precision spectra of the primary 
component. 



5 BINARITY AND EVOLUTIONARY STATE 
5.1 Binarity 

Determining whether the HD 72106 system is truly a binary 
or just an optical double star (i.e. an accidental conjunction 
of 2 stars at different distances along the line of sight) is 
critical. If the system is a coeval binary, this fact allows us 
to constrain the age of the primary much more accurately 
than would otherwise be possible. Additionally, it suggests 
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Figure 3. Position angle measurements of the HD 72106, ex- 
tracted from the WDS, together with their best fit line. 



that both components formed from approximately the same 
material, making HD 72106 an interesting system from the 
point of view of stellar magnetic and chemical evolution. 

The system has a projected separation of 0.805 arcsec, 
and is thus fairly widely separated. Given the Hipparcos 
parallax of 3.60± 1.14 mas, this implies a minimum physical 
separation of 224 ± 71 AU. 

The Hipparcos observations were solved as a binary 
system producing a "good quality" solution for the sys- 
tem. This implies that the stars have the same parallax 
and proper motions, at the precision of the Hipparcos ob- 
servations. Hipparcos finds a large proper motion in right 
ascension of —6.41 ± 0.82 masyr" 1 and in declination of 
7.97± 1.23 masyr -1 l|van Leeuwenll2007al ). In the spectrum 
fitting procedure, described in Section [7] we included he- 
liocentric radial velocities for each star as a free parameter. 
From this fitting, we found identical radial velocities of 22±1 
kms -1 for both components. Thus it appears that the stars 
are at the same point in space and moving together in three 
dimensions, strongly suggesting that they are in fact physi- 
cally associated. 

Dr. Brian Mason at the United States Naval Obser- 
vatory kindly provided a record of separation and position 
angle observations from the Washington Double Star Cat- 
alog (WDS), dating back to 1902 (private correspondence, 
Mason, 2008). These have been compiled from a number 
of sources in the literature, some of which are no longer 
readily accessible. These observations indicate no significant 
change in the separation of the components in 89.37 years, 
but show a clear and systematic increase in the position 
angle of the stars of ~10.5°, illus trated in Fi g ure [3] This 
relative motion was not noted by IWade et al l (|2007r ). pre- 
sumably because they examined only a subset of these ob- 
servations. The change in position angle occurs at a rate of 
0.117±0.007° yr _1 . In contrast, there is no trend in the set 
of separation measurements, with an average separation of 
0.786 arcsec and a standard deviation of 0.074 arcsec. This 
kinematical behavior, together with the identical radial ve- 
locities and consistent parallaxes of the components, makes 
it high likely that the system is a true binary. 

For the system to be a true binary it must be gravita- 
tionally bound. We have therefore verified that the observed 
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properties are consistent with binarity. We assumed circu- 
lar orbits, that the minimum possible separation is the true 
separation (224 ±71 AU), and used the evolutionary masses 
derived in Section T5.2I These assumptions are necessary, as 
we cannot determine the true separation and orbit of the bi- 
nary without much more accurate observations over a much 
longer period of time. In this geometry, the relative veloc- 
ity of the stars is 4.1 ± 0.7 kms -1 and the orbital period 
is 1600 ± 600 years. This relative velocity is larger than ob- 
served, implying that the orbital plane is not parallel to 
the line of sight, or that the stars' physical separation is 
larger than we have assumed (which it almost certainly is). 
If the orbital plane of the stars is in the plane of the sky, 
the orbital period implies that an observed change in posi- 
tion angle of 0.219 ± 0.079° yr _1 should be seen, which is 
roughly consistent with the WDS observed rate of change. 
Thus we conclude that a wide range of orbits are consis- 
tent with the observed radial velocities and position angle 
change. One example of an orbit that is consistent with all 
observations is a circular orbit with a semimajor axis of 350 
AU (implying a period of 3200 years), with an inclination 
of the orbital plane to the line of sight of 44.5°, and a po- 
sition angle of the orbital axis of 200°. In this scenario, the 
two components would both be crossing the line-of-sight at 
the current epoch, yielding zero relative radial velocity. The 
change in position angle would be 0.112° yr _1 and a ap- 
parent separation would be 0.807 arcsec. In this example, 
in 800 years the components would achieve their maximum 
difference in relative radial velocity of 2.3 km/s, and their 
maximum apparent separation of 1.151 arcsec. Of course, 
there is no reason to prefer this particular solution over the 
many others that are consistent with observations. 

When placed on the Hertzsprung-Russell (H-R) dia- 
gram, as discussed in Section f5. 2 1 the stars are found to have 
positions that allow for a range of coeval solutions. Thus ac- 
cording to the H-R diagram positions, it is quite possible 
for the stars to have formed and evolved together, further 
supporting the binary hypothesis. 

Thus we conclude that HD 72106 is very likely a true 
binary system. The stars are at the same position in space, 
moving in the same (three dimensional) direction at the 
same velocity, allow for a wide range of gravitationally 
bound orbits, and appear to be consistent with coeval evo- 
lution. 

5.2 Fundamental Physical Properties 

Effective temperature and surface gravity for both com- 
ponents of HD 72106 were determined by fitting Balmer 
lines. Unfortunately, no intermediate-band photometry of 
the components of HD 72106 was available. The Ha Balmer 
line in HD 72106B is clearly contaminated with emission, 
and hence was not considered in this analysis. The other 
Balmer lines of both stars appear to be free of emission, 
and lines of each star can be well-fit with a single model 

atm osphere. 

I Wade et al.l (|2005l ) performed Balmer line fitting for 
both components of HD 72106 using observations from the 
FORS1. For the primary they found: T cff = 11000 ± 1000 
K and 3.5 < log g < 4.5, (best fit T cff = 11000 K log g = 
4.0). For the secondary they found T cff = 8000 ± 500 and 
4.0 log g sS 4.5, (best fit T cff = 8000 K log g = 4.5). 



We repeated the fitting procedure with the FORS1 spectra, 
using the method outlined abo ve. Solar abun dance ATLAS9 
model atmospheres were used (|Kuruc3[l993) to produce the 
synthetic Balmer lines and the model lines were convolved 
with a Gaussian instrumental profile of appropriate width 
to match the observ ations. We arrived at results identical 
to I Wade et all i|2005l ) for the primary. For the secondary we 
find the temperature range 7500 K to 9000 K and the log g 
range 4.0 to 4.5 more realistic, with a best fit value of 8000 
K at log g— 4.5. These values from the FORS1 observations 
are also consistent with those determined from Balmer line 
fitting of our ESPaDOnS spectra. 

In Section [7] we performed detailed spectrum synthesis 
of both components of HD 72106 and found that the atmo- 
spheric parameters of the secondary derived from Balmer 
lines are not compatible with the metallic line spectrum. In 
particular, we were unable to simultaneously fit lines of Fe I 
and Fe II, Cr I and Cr II, and Ti I and Ti II, suggesting a 
problem with the local thermodynamic equilibrium (LTE) 
ionization balance of these species. Including T e ff and log g 
as free parameters in the spectrum fit allowed us to sat- 
isfactorily match the metallic line spectrum, provided that 
Tcff ~8750 K and log g ~4.0. While these values produce 
a somewhat poorer fit to the Balmer lines, the fit to metal- 
lic lines is improved substantially. Hence we conclude that 
the higher temperature and lower log g provide a better de- 
scription of the atmosphere of HD 72106B. 

Including the full range of T c g and log g which pro- 
vide acceptable fits to both the metallic line spectra and the 
Balmer lines of HD 72106B, we adopt T e f£ = 8750 ± 500 K 
and log g = 4.0 ± 0.5. 

Using the Hipparcos para llax of the HD 72106 sys - 
tem and the Tycho magnitudes (jFabricius fe Makarovll200oh 
we determined the luminosities of both components of HD 
72106. To convert the observed Tycho V magnitudes (Vt = 
9.00 ± 0.01 for the primary and V T = 9.62 ± 0.02 for the 
secondary) into a Johnson V magnitudes (V) we used the 
empirical relation \ES All 19971 ): 

V = V T + 0m(B T - V T ), (3) 

where Bt is the Tycho B magnitude. For the secondary, 
the bolomet ric correction for main sequence stars from 
iGravl (120051 ) was used, yielding a value of 0.01 ± 0.06 
m ag. For the primary, th e bolometric correction relation 
of lLandstreet et all (|2007l ) was used, yielding a value of 
—0.37 ± 1.9 mag. This calibration is tailored specifically 
for magnetic chemically peculiar A and B type (Ap and 
Bp) stars. As will be shown, the primary has properties 
ve ry similar to those of a Bp star, thus the calibration 
of lLandstreet et al.l l|2007l ) is more appropriate. From this 
procedure we find that the luminosity of the primary is 
22 Z/q, and the luminosity of the secondary is 9.2 +|'g Lq. 
From the Stefan-Boltzmann equation we find identical radii, 
given our uncertainty, of 1.3 ± 0.5 -Rq for the primary and 
1.3 ± 0.5 Rq for the secondary (with la error bars). 

With the luminosity and effective temperature deter- 
mined, we can place the stars on the H-R diagram, as shown 
in Figure [4] and compare their positions with theoretical 
evolutionary tracks and isochrones. Using the pre-main se- 
quence e volutionary m odel calculations fro m CESAM (ver- 
sion 2 K; iMorell 119971 1 and the birthline of iPalla fe Stahlerl 
(|l993l ). we find the mass of the primary to be 2.4 ± 0.3 Mq 
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Table 2. Fundamental physical properties of HD 72106A and B. 
Magnetic field parameters are determined in Section [8] 



HD 72106A HD 72106B 



T cff (K) 11000 ±1000 8750 ±500 

log g (cgs) 4.0 ±0.5 4.0 ±0.5 

L(L Q ) 22120 9.2+8.2 

R(Rq) 1.3 ±0.5 1.3 ±0.5 

M (Mq) 2.4 ±0.3 1.9 ±0.2 

age (Myr) 6-13 6-13 

usini (kms" 1 ) 41.0 ±0.3 53.9 ± 1.0 

P rot (days) 0.63995 ± 0.00009 

i(°) 24 ±10 

/3(°) 57 ±5 

Bp (G) 1230 ± 80 



Figure 4. An H-R diagram containing both components of the 
HD 72106 system. Isochrones (solid lines), evolutionary tracks 
(dashed lines), and the zero age main sequen ce (ZAMS) line, are 
calculated from CESAM models dMorelll 19971) . T he birth line (for 
an acc retion rate of 10 — 5 Mq yr — 1 ) is taken from lPalla &: Stahlerl 
dl993l l. 



and mass of the secondary to be 1.9±0.2 Mq (with la error 
bars). From this we derive an evolutionary log g of 4.6 ±0.3 
for the primary, and 4.5 ± 0.3 for the secondary, both of 
which are consistent with (although systematically higher 
than) the spectroscopic values. We determine the binary's 
age to be between 6 and 13 Myr (measured from the birth 
line), based on the position of the secondary in the H-R 
diagram and assuming a protostellar accretion rate of 10 -5 
Mq yr" 1 . In doing this we have assumed, as most investiga- 
tors do, that the presence of circumstellar material implies 
that the HAeBe secondary is on the pre- main sequence, and 
thus that we can constrain its age with pre-main sequence 
evolutionary tracks. Given the derived admissible range of 
age and mass, the primary could still be on the pre-main se- 
quence, reaching the zero-age main sequence (ZAMS) in less 
than 1 Myr. The most likely case (from 'best fit' positions) 
is that the system is ~10 Myr old, and the primary has just 
entered the main sequence, while the secondary is on its final 
approach to the ZAMS. In this case the primary would have 
spent ~6 Myr on the main sequence, giving it a fractional 
age on the main sequence (r) of 0.01. In the oldest limiting 
case, the secondary is just reaching the ZAMS and the pri- 
mary has been on the main sequence for ~9 Myr giving it a 
fractional age of 0.015. Further observations, particularly a 
more accurate distance measurement, are necessary to more 
precisely determine the evolutionary status of the primary. 
Regardless, while HD 72106A may not be on the pre-main 
sequence, it is certainly one of the youngest known main 
sequence stars of its type. 

The physical properties we derive for both stars are 
summarized in Table [2] It is worth noting that, while the 
absolute luminosities of the components are poorly deter- 
mined, their ratio (La/ Lb = 2.3 ± 0.4) is very well deter- 
mined. This is because the major uncertainty in the absolute 
luminosity is the distance to the stars, and the stars are lo- 
cated at the same distance. Thus the values of the radii or 
masses derived for the two stars are not independent, since 



spacing between the components in logL on the H-R dia- 
gram must remain fixed. 



6 LEAST SQUARES DECONVOLUTION AND 
STELLAR ROTATION 

6.1 LSD and Longitudinal Magnetic Field 

In order to measure the magnetic fields of HD 72106A 
fc B, we employe d Least Squares Deconvolution (LSD; 
iDonati et al.l 11997m . This cross correlation technique uses a 
table of input atomic data (the line mask) to produce a 
deconvolved mean line profile from an observed spectrum. 
Using a large number of lines, ~7000 for our spectra, can re- 
sult in dramatic improvements in S/N. LSD was performed 
on all our observations of the individual components of HD 
72106, as well as the reconstructed spectra of HD 72106A. 
The full set of LSD profiles in Stokes I and V for HD 72106A, 
calculated using specific elements, are presented in Section 
[S]in Figure [TO] 

The line mask used for the primary was calculated as- 
suming an effective temperature of 11000 K, log g of 4.0, 
typical Ap star abundances, and a line depth cutoff of 0.1 
(as a fraction of the continuum). This produced LSD profiles 
with a mean wavelength of 523.4 nm, a mean Lande factor of 
1.27, and a mean excitation potential of 5.97 eV, with small 
variations between observations. These mean values repre- 
sent the average values of the lines used in the computation 
of the LSD profile, weighted by the S/N of the lines. The 
computed LSD profiles are relatively insensitive to variations 
or errors in the line mask in temperature or abundance (e.g 
IShorlin et al1l2002r i. For the secondary, a T eS = 8750 K, 
log g = 4.0 line mask was used with solar abundances and 
a line depth cutoff of 0.1. With this mask we found a mean 
wavelength of 533 nm, a mean Lande factor of 1.21, and a 
mean excitation potential of 3.99 eV. 

The longitudinal magnetic field was then measur e d from 
each set of LSD profiles, using Eq. (1) of lWade et aD l|2000h . 
Integration was performed through the portion of the line 
profile that exceeded 15% of the total line depth, with an ad- 
ditional 5 kms -1 on either end. The longitudinal field data 
are presented in Table Q] We find a maximum longitudinal 
field for HD 72106A of 374 ± 32 G. 
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In the secondary, in contrast, we detect no longitudinal 
field, with error bars of 50-150 G. Moreover, the LSD profiles 
of HD 72106B yield no e vidence of a magne tic field. The 
analysis of HD 72106B by IWade et~ai] (|2007h . discussed in 
Section [2] found no magnetic field in either the Balmer lines 
or metallic lines with error bars of ~100 G Thus we conclude 
that if the secondary has a magnetic field, its longitudinal 
component never exceeds ~ 200 G. 

6.2 Rotation Period of HD 72106A 

Variability i n the spectrum of HD 72106A was noted by 
IWade et all (|2005T ). A careful reexamination of their data 
shows that most of their reported spectrum variability was 
due to the lack of atmospheric dispersion correction during 
the acquisition of one of their spectra. However, we do see 
substantial metal line variability in all of our spectra of the 
primary, with a variation timescale on the order of a day. 
We do not confir m the rotational pe riod of approximately 2 
days proposed by I Wade et ail (|2005T ). 

In order to investigate whether this variability in HD 
72106A could be due to rotational modulation, we fit the 
longitudinal magnetic field observations with a first order si- 
nusoid of fixed period, but free phase, amplitude, and mean. 
This process was repeated for many periods to construct 
a period ogram in a fashion analogous to the Lomb-Scargle 
method l|Press et all 1 19921 ). In the resulting periodogram, 
shown in Figure [6] (left frame) , we found a number of peri- 
ods that provided a good fit to our observations, with the 
best period at about 0.640 days. Our longitudinal field ob- 
servations, phased with this best fit period, are shown in 
Figure [5] We interpret this variability as being due to the 
rotation of HD 721 06A, in the context of the oblique rotator 
model (ORM; e. dStibbsll95Cl ). Using the radius determined 
in Section r5.2l together with the best fit period, we find that 
a range of inclination angles will satisfy the derived v sin i 
from Section r7.ll The smooth sinusoidal behavior of the lon- 
gitudinal magnetic field, the stability of the period over two 
years, and the period consistent with the star's v sin i and 
radius all support our interpretation of the variability in the 
context of the ORM. It is worth noting that the phase cov- 
erage for periods around 0.64 days is poor between phases 
0.4 and 0.7. However, the remainder of our analysis is based 
on rotationally broadened LSD profiles rather than longitu- 
dinal field measurements, which mitigates this phase gap. 

In order to search further for periodic variations in the 
observations of HD 72106A, we examined the variability 
within the Stokes / and V LSD profiles. To do this, we 
constructed a periodogram based on the variability of each 
pixel of the LSD profiles, using the same technique as for the 
longitudinal field measurements. The periodograms for each 
pixel were then averaged and weighted by the amplitude of 
variability in that pixel, to produce a periodogram for the 
entire set of LSD profiles. This procedure was performed for 
both Stokes I and Stokes V LSD profiles, and the resulting 
periodograms are shown in Figure [6] (right frame) . 

Periods at 0.38983 days, 0.63995 days, 1.6921 days and 
1.7859 days all produced minima in one or more of the pe- 
riodograms that were within the 99% confidence limit of 
the global minimum. These candidate periods were each in- 
vestigated by eye. Each period was used to determine the 
phasing of the LSD profiles, and the profiles were exam- 




Figure 5. Longitudinal magnetic field measurements of HD 
72106A (points), phased with the adopted 0.63995 day period, 
and the longitudinal field curve from the best fit dipolar mag- 
netic geometry (solid line) determined in Section [8] 



ined for variability physically consistent with rotation. For a 
physically sensible rotational variation, features must move 
smoothly across the profile from blue to red, and there must 
not be large changes in the profile with small changes in 
phase. Upon careful examination, the period of 0.63995 days 
was the only period to produce a physically sensible phas- 
ing of the I and V profiles. Thus we adopt the rotational 
ephemeris (with maximum longitudinal field at zero phase) 
for HD 72106A: 

HJD = (2453747.017 ± 0.013) + (0.63995 ± 0.00009) ■ E. (4) 

The LSD profiles phased according to this ephemeris are 
shown in Figure QIJ] This period is notable in that it is 
one of the shortest rotation periods seen in any magnetic 
intermediate-mass star. 



7 SPECTRUM SYNTHESIS 

The ZEEMAN2 spectrum synthesis code was used to model 
the observed spectra of HD 72106A & B. This code solves 
the polarized radiative transfer equations in local thermo- 
dynamic equilibrium. In order to find the best fit model 
of the observed spectrum, a Levenberg-Marquardt \ 2 mini- 
mization routine was implemented. This routine consistently 
produced high quality, stable solutions with good efficiency. 
The results were, in all cases, carefully examined by eye to 
verify that they represented both a physical solution and the 
true best fit, rather than a local minimum in \ 2 ■ 

Chemical abundances (in the form [^y 2 -] , assuming 
Nne/Nu = 0.098), projected rotational velocity (usini), 
and for HD 72106B, microturbulence, were the free param- 
eters in the fit. Chemical abundances were assumed to be 
uniform both vertically and horizontally. This is not strictly 
true for HD 72106A, which displays patchy distributions 
of some elements, discussed further in Section [8] Thus the 
abundances derived for HD 72 106 A are averages across the 
visible disk of the star at the phase of the observation used 
(the 12 Jan. 06 observation, rotation phase 0.53). However, 
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Figure 6. Periodograms for HD 72106A based on longitudinal field measurements (left panel) and LSD profile variability (right panel) 
in Stokes / and V. A few minima in each periodogram fall below the 99% confidence limit (horizontal line), however the deepest minimum 
is at 0.64 days in all cases. 



since there is not much change in equivalent widths be- 
tween spectra, the observation of HD 72106A used is approx- 
imately representative of the global average for all elements, 
even those with patchy distributions. 

Microturbulence is not expected in the primary, due to 
its strong magnetic field. By analogy to Ap/Bp stars, the 
magnetic field would likely suppress microturbulence, thus 
microturbulence was fixed at kms" 1 for HD 72106A. This 
did not produce any detectable discrepancies between our 
model and observed spectra. For the primary we adopted the 
magnetic field model derived from the Stokes V LSD pro- 
files in Section [S] Uncertainties in this magnetic field model 
have a negligible impact on the derived abundances. For the 
secondary no magnetic field was included. 

Standard model atmos pheres with s olar abundance, 
computed with the ATLAS9 l)Kurucdl 19931 ) code, were used 
as input for ZEEMAN2. Input atomic line data was ob- 
tained from the V ienna Atomic Line Database (VALD) 
(|Kupka et all I1999T ). using the "extract stellar" utility in 
VALD (with the default line list configuration). 

Seven segments of spectrum, ranging from 100 to 200 A 
in length, were independently fit for each star (4170-4265 A, 
4400-4500 A, 4500-4600 A, 4600-4700 A, 5000-5200 A, 5200- 
5400 A, and 5400-5600 A) The final best fit parameters were 
taken as the average of the parameters derived from indi- 
vidual windows. Uncertainties were taken as the standard 
deviation of the parameters. For abundances that were only 
constrained in a few windows, the uncertainties were esti- 
mated by eye, by taking the change in abundance necessary 
to shift the synthetic spectrum beyond the noise in the ob- 
served spectrum and any possible continuum normalization 
errors. In these case the presented uncertainties are indi- 
cated by an asterisk (*) in Table [3] 



7.1 Abundances in HD 72106A 

Sample fits of our synthetic spectra to the observations of 
HD 72106A are presented in Figure [7] The average best 
fit abundances and vsini, with uncertainties, are given in 
Table [3] The abundances are presented graphically, relative 
to solar abundances, in Figure [5] 

Remarkably strong overabundances of Si, Cr, Fe and Nd 
are found in HD 72106A. The Si and Fe abundances appear 
to be above solar by ~1 dex, and Cr appears to be enhanced 
by ~2 dex, whereas He appears to be ~1.5 dex underabun- 
dant. A number of elements, such as Al, Sc and Sr appear 
to have solar abundances. A couple of elements, particularly 
Mg and O, hint at possible peculiarity but require further 
study before concrete conclusions can be drawn. We find 
vsini = 41.0 ± 0.7 kms -1 , which is fairly low for a main 
sequence B star, but within the normal range for Bp stars. 
The strong overabundances in Si, Cr, Fe, and Nd, as well as 
the unde rabundance in He, are co mmon features of cooler 
Bp stars (jjaschek fc Jascheklll995h . 

7.2 Abundances in HD 72106B 

HD 72106B was initially modeled with an effective temper- 
ature of 8000 K and log g = 4.5. However, as discussed in 
Section [5.21 this model provided a poor fit to the observed 
stellar spectrum, and hence T c g and log g were included 
as free parameters in our fit. Chemical abundances, vsini 
and microturbulence were also included as free parameters 
in the fit, while the magnetic field was set to zero. From this 
procedure we found a best fit T c s of 8750 ± 500 K and log g 
— 4.0 ± 0.5. Sample best fit synthetic spectra can be seen in 
Figure [7J compared with the observed spectrum used in the 
fitting process. The average best fit abundances, vsini, and 
microturbulence derived for HD 72106B are presented in Ta- 
ble [3] and the abundances are shown graphically relative to 
solar abundances in Figure [8] 



10 Folsom et al. 




>-o CZ i I i I i I i I i ZJ I i I i I I I i I i I 

5000 5010 5020 5030 5040 5050 4200 4210 4220 4230 4240 4250 

Wavelength (angstroms) Wavelength (angstroms) 



Figure 7. Sample best fit synthetic spectra for HD 72106A (left panels) and HD 72106B (right panels) in two independently fit spectral 
windows. The observations of the individual stellar components from Mar. 5, 2007 are used for both stars. Major contributors to each 
line have been labeled, in order of importance. The smooth solid line is the best fit spectrum in this region, the dashed line is a spectrum 
computed with solar chemical abundances. Dots represent the observations. 



The large majority of elements are consistent with so- 
lar abundances, within 2a at most. A few elements appear 
to depart marginally from solar values, with a significance 
slightly greater than 2a. C appears to be overabundant with 
~ 2a significance, whereas Sc is underabundant by ~ 3a. 
Under abundances of Sc and Ca are characteristic of Am 
stars. However, since the Ca abundance is approximately 
normal in HD 72106B, and iron peak elements are not over- 
abundant, HD 72106B is not Am star. Thus, nearly all the 
elements are within 2a of solar abundance, and no elements 
display the strong peculiarities seen in HD 72106A. This re- 
sult is consistent with the approximately solar abu ndances 
seen in most HAeBe stars |Acke fe Wae lkcns 2004). It is of 
possible importance that the best fit abundances are con- 
sistently below solar, albeit with small significance levels. 
However there remains a large uncertainty in our best fit 
temperature and log g. An increase in the adopted temper- 
ature by 250 K to 9000 K would systematically increase the 
best fit abundances by 0.1 dex or 0.2 dex, essentially elimi- 
nating this trend. 



8 MAGNETIC DOPPLER IMAGING 

The observed metallic line variability with rotation strongly 
suggests that there are horizontal chemical abundance in- 
homogeneities in the atmosphere of HD 72106A. In an ef- 
fort to determine the structure of these inhomogeneities, as 
well as the magnetic field structure, we performed Magnetic 
Doppler Imaging (MDI) on HD 72106A. Doppler Imaging is 
a method for inverting a timeseries of variable line profiles 
at known rotational phases in order to reconstruct the sur- 
face distributions that give rise to the observed variations. 
In MDI, polarized line profiles are also included in the in- 
version, allowing for the reconstruction of a magnetic field 
through the Zeeman effect. 

We used the MDI code INVERS10, developed by O 
Kochukhov and N. Piskunov (jPiskunov fc Kochukhovll2002l ; 
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Figure 8. Abundances relative to solar for HD 72106A (circles) 
and HD 72106B (squares), averaged over all spectral windows 
m odeled. The dashed lin e at represents solar abundance, based 
on lGrevesse et al. Points marked with an arrow indicate 

the value is an upper limit only. Strong departures from solar 
abundance can be seen for HD 72106A, whereas HD 72106B has 
largely solar abundances. 

I Kochukhov fc Piskunovl 120021 ) for surface chemical abun- 
dance mapping and magnetic field reconstruction. This 
program performs accurate polarized LTE spectrum syn- 
thesis, using pre-calculated model stellar atmospheres. IN- 
VERS10 allows for simultaneous modeling of multiple chem- 
ical elements and multiple wavelength regions, and takes 
into account blended lines. Tikhonov regulari zation is used 
l|Tikhonovlll963l : IPiskunov fc Kochukhovl2002l ) , to constrain 
abundance gradients in the surface map. While Tikhonov 
regularization can be applied to the mag netic field as well, 
we o pted for multipolar regularization (|Kochukhov et aD 
2002), due to the lack of Stokes Q and U spectra and the 
relatively low S/N in the Stokes V spectra. 
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Table 3. Averaged best fit chemical abundances, ?;sini and mi- 
croturbulenc e (g) for HD 7 2106A and B as well as solar abun- 
dances from iGrevesse et al, | j2005h . Entries marked by an * are 
based on only a few lines. 





HD 72106A 


HD 72106B 


Solar 


v sin i (km s 1 ) 


41.0 ± 0.6 


53.9 ± 1.0 




£ (kms- 1 ) 




2.3 ± 0.6 




He 


-2.8 i 0.2 * 




-1.07 


Q 






-3.61 





-3.0 ± 0.15* 




-3.34 


IVlg 




a fin -t- n i (\ 


A A7 

-41. "4: ( 


Al 


-5.8 ± 0.2 * 




-5.63 


Si 


-3.60 ± 0.14 


-4.7 ± 0.3 


-4.49 


C'a 


-5.1 ± 0.3 * 


-6.0 ± 0.2 


-5.69 


Sc 


-8.6 ± 0.2 * 


-9.13 ± 0.08 


-8.83 


Ti 


-6.04 ± 0.10 


-7.24 ± 0.09 


-7.10 


Cr 


-4.33 ± 0.13 


-6.32 ± 0.19 


-6.36 


Mn 




-6.7 ± 0.3 * 


-6.61 


Fe 


-3.49 ± 0.07 


-4.64 ± 0.17 


-4.55 


Ni 




-6.3 ± 0.3 


-5.77 


Sr 


-8.7 ± 0.4 * 


-9.2 ± 0.4 * 


-9.08 


Y 




-10.0 ± 0.3 


-9.79 


Ba 


< -8.5 * 


-10.2 ± 0.4 * 


-9.83 


Ce 




< -9.0 * 


-10.30 


Nd 


-7.5 ± 0.3 * 


< -9.2 * 


-10.54 



Input into the MDI model of HD 72106A included the 
adopted effective temperature and surface gravity from Sec- 
tion [ITSI T cft = 11000 K and log g = 4.0. The inclina- 
tion angle of the star's rotation axis to the line of sight (i) 
was calculated using the radius determined in Section 15.21 
the rotation period from Section 16.21 and the v sin i deter- 
mined from spectrum fitting in Section [7. II giving an angle 
of i = 24 ± 10°. The phasing of observations was deter- 
mined from the adopted ephemeris. An initial vsini of 41 
kms -1 and abundances, for the treatment of blended lines, 
were taken from the results of our abundance analysis of HD 
72106A, presented in Section \7. II The initial magnetic field 
geometry was assumed to be a dipole. 

Initially, five Si II lines at 4128 A, 4130 A, 5056 A, 5041 
A, and 6371 A were considered for Doppler Imaging. How- 
ever, the S/N of the observations was insufficient to produce 
high quality Doppler maps for Si, as well as for all other el- 
ements. 

In order to improve the S/N of our observations we cal- 
culated LSD profiles for individual chemical elements, and 
then performed MDI using those profiles. The use of LSD 
profiles to determine magnetic field geometries, stellar p ul- 
sation, and surface structures (|Donati et al 1 l2000l .l 20011 ) is 
well established. LSD profiles of individual elements were 
obtained by creating LSD line masks containing only lines 
of the element of interest. The line masks used were derived 
from the 11000 K Ap star line mask discussed in Section 
[U Doppler Imaging was then performed using the set of 
LSD profiles for each element as described above. The el- 
ements Si, Ti, Cr, and Fe were used in this process. Due 
to the low S/N in Stokes V for the purposes of MDI, even 
in the LSD profiles, the magnetic field geometry was de- 
termined from the higher S/N Cr and Fe profiles, and held 



Table 4. Mean atomic data for the element specific LSD profiles 
used in Magnetic Doppler Imaging. 





Si 


Ti 


Cr 


Fe 


mean wavelength (A) 


5580 


4970 


5190 


5310 


mean excitation (eV) 


7.68 


2.03 


5.65 


6.46 


log mean gf 


-0.17 


-0.63 


-0.26 


-0.35 


mean Lande factor 


1.35 


1.12 


1.29 


1.27 



fixed for the lower S/N Si and Ti profiles. Additionally, the 
multipolar magnetic field regularization was restricted to to 
I = 1 modes, effectively providing a 'dipolar regularization'. 
Mean LSD profile atomic data were used for the Doppler 
Imaging process, providing wavelengths and excitation po- 
tentials, shown in Tabled Since the depth of the LSD pro- 
file is a complex function of the lines used in the analysis, 
the mean gf value may not represent a realistic oscillator 
strength for the LSD profile. As a consequence, the absolute 
abundance scale of the Doppler images is somewhat uncer- 
tain. 

The final abundance maps, and magnetic modulus and 
vector maps, are shown in Figure [9] The corresponding best 
fits to the LSD line profiles are presented in Figure 1101 In 
this process we have assumed that the LSD profile behaves 
like a real spectral line, and that the mean atomic data is 
approximately representative of the LSD profile. Addition- 
ally, unaccounted for line blends in the LSD process could 
slightly distort the LSD profiles, adding some uncertainty to 
the finer details of the Doppler map. Thus, while the large 
scale features of the Doppler Imaging process are fairly cer- 
tain, one should use some caution in interpreting the finer 
details of the map, as well as the absolute abundance scale. 
Similarly, the use of LSD profiles to reconstruct the mag- 
netic field geometry through MDI could potentially intro- 
duce some systematic uncertainty into the derived geome- 
try. However, the magnetic geometry derived through MDI 
is fully consistent with the longitudinal magnetic field vari- 
ability and nearly identical to that we derive by directly 
modeling the LSD Stokes V profiles with a simple dipole. 

Strong inhomogeneities are reconstructed for all four el- 
ements. Ti, Cr, and Fe all seem to share similar abundance 
patterns, though Cr appears to have somewhat larger, more 
elongates spots. The similar distributions are reflected in the 
similar phase variations of the LSD profiles in Figure 1101 A 
large patch of overabundance centered near phase is appar- 
ent in all three maps, with another somewhat smaller over- 
abundance spot about 180° away in longitude, at the same 
latitude, around phase 0.6. The Si map shares the larger 
spot but not the smaller, and this is reflected in the Si LSD 
profiles. There appears to be a large overabundance spot of 
Si at the equator (seen near the limb) near phase 0, although 
the sensitivity of the map is poor in that region due to its 
small projection. Note that the fits to the relatively noisy 
Si LSD profiles are somewhat poorer than the fits for the 
other elements analyzed, thus the results for Si are some- 
what more uncertain. The magnetic field geometry derived, 
shown in Figure [9] is essentially described by a dipole with 
a strength at the magnetic pole (-Bp) of 1230 ± 80 G, and 
an obliquity angle of the magnetic dipole with respect to 
the rotation axis (/?) of 57° ± 5. This dipole magnetic field 
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Figure 9. Surface maps of Si, Ti, Cr, and Fe abundances, as well as the magnetic field, for HD 72106A. The maps are all based on fits 
to LSD profiles. The 'X' represents the rotational pole, the thick line circle indicates the rotational equator. The scale on the right in 
the abundance maps is in units of log jj-^-- The map of the magnetic field intensity is labeled 'a' and the magnetic field direction is 'b' 
(shown as arrows). 



geometry is in full agreement with the observed longitudinal 
magnetic field variability. If there are any departures from 
a purely dipolar field geometry, they are not evident in our 
data. When the magnetic field geometry is compared to the 
abundance maps, it appears that the positive magnetic pole 
lies near the large spot of overabundance at phase in all 
four maps. However, the magnetic pole is offset from the 
center of the abundance spot, being nearer the rotational 
equator and at a slightly different latitude, thus the rela- 
tionship is not entirely clear. 

Interestingly, the pattern of abundances seen here for 
HD 72106A bear some similarities to th ose obtained for e 
Ursae Majoris by iLueftinger et al.l (|2003l ). e Ursae Majoris 
(HD 112185) is a 9000 K (log g = 3.6) late main sequence Ap 
star with ~ 5 d ay period and a dipole fi eld strength of several 
hundred gauss. ILueftinger et all (|2003T ) constructed Doppler 
maps of Ti, Cr and Fe, as well as Ca, Mg, Sr, and Mn. They 
found distributions of Cr and Fe very similar to each other, 
with two large spots of overabundance near the longitude of 



the magnetic poles. Ti was roughly anti-correlated with Fe 
and Cr, displaying two large spots of underabundance at the 
same positions as the overabundance spots of Cr and Fe. We 
see little variability in Ca, and have insufficient S/N in Ca, 
Mg, Sr, and Mn to construc t Doppler maps for comparison 
with ILueftinger et all (|2003h . The results for Fe and Cr are 
qualitatively very similar to what we see for HD 72106A, 
however the results for Ti differ significantly. The signifi- 
cance of these similarities and differences is unclear, as the 
evolution of chemical abundance spots over a star's main 
sequence lifetime is not understood. 



9 DISCUSSION AND CONCLUSIONS 

We have analyzed 20 high resolution spectropolarimetric ob- 
servations of the HD 72106 system. In these observations we 
see clear evidence for a magnetic field in HD 72106A. We also 
confirm that HD 72106B is a HAeBe star, based on emission 
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Figure 10. Fits of synthetic Stokes I (for Si, Ti, Cr, and Fe) and Stokes V (for Cr and Fe) LSD profiles to observed LSD profiles. The 
profiles are labeled according to element, phase and mean wavelength. The bars near the bottom of the figure indicate the vertical and 
horizontal scale, 0.5% of the continuum and 0.5 A respectively. Generally good fits can be seen, ranging from the high S/N Fe profiles 
at the best to the noisy Si profiles at the worst. 



in Ha and the O I 7771 A triplet, and that it displays no 
magnetic field. 

There is strong evidence that the HD 72106 system is 
a true binary system. Both stars have the same proper mo- 
tions, the same radial velocities, and Hipparcos solution to 
the system places both stars at the same distance. Although 
the separation of the components has remained constant in 
the pas 90 years, there has been a slow systematic increase 
in the position angle observed. Additionally there is a wide 
range of possible bound orbits consistent with all available 
observations. 

We find the age of the system to be between 6 and 13 
Myr, based on the H-R diagram position of the secondary 
(and assuming that, as a HAeBe star, the secondary is a 
pre-main sequence star). Thus the system is fairly evolved, 
for a pre-main sequence system. In the youngest limit, the 



primary would within 1 Myr of the ZAMS. However, in the 
oldest limit the primary would have reached the main se- 
quence 9 Myr ago, but it would have only passed about 
1.5% of its main sequence lifetime. Thus, while it may not 
be on the pre-main sequence, HD 72106A is certainly very 
young. While it is possible that the system is not coeval and 
instead was form by capture, the data available (consistent 
H-R diagram positions and the young age of the secondary) 
suggests that the system truly coeval. Even with the large 
uncertainty in age, only a few known Ap/Bp stars approach 
the maximum fractional age of HD 72 106 A l|Bagnulo et al.l 
|2004| ; iLandstreet et alj|2007l ). Given its evolutionary status, 
HD 72106A appears to represent a link between magnetic 
HAeBe stars and the Ap/Bp stars. 

HD 72106A possesses a strong, predominantly dipolar, 
magnetic field. We find that a centered dipole with a po- 
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lar field strength of 1230 ± 80 G and an obliquity angle 
of 57 ± 5° is sufficient to provide a good fit to our observa- 
tions. In HD 72106B we find no evidence of a magnetic field, 
with an upper limit on the longitudinal field of about 200 
G. HD 72106A is one of the very youngest stars for which 
detailed modeling of the magnetic field geometry has been 
performed. The only clea rly younger ob j ect is t he Herbig Be 
star HD 200775A, which lAlecian et alj (|2008ah find to have 
an age of 0.1 ± 0.05 Myr. This star also displays a dipolar 
magnetic field geometry, with a magnetic field strength at 
the pole of 1000 ± 150 G and an obliquity of 12 5 ±8°, slightly 
offset at 0.05 ± 0.04 _R* from the star's center (|Alecian et al.l 
1200881 ). HD 200775A and HD 72106A have similar magnetic 
field characteristics, both of which are analogous to most 
main sequence Ap/Bp stars. This provides further evidence 
that there is a continuum of magnetic A and B stars from 
the pre-main sequence through the main sequence. 

We see strong chemical peculiarities in HD 72 106 A, par- 
ticularly an underabundance of He and overabundances of 
Si, Ti, Cr, Fe, and Nd. Thus, from the point of view of abun- 
dance analysis, HD 72106A appears to be a Bp star. This 
implies that chemical peculiarities can form very early in a 
star's lifetime, possibly even on the pre-main s equence. This 
conclus ion is supported by the recent results of lAlecian et al.l 
(2008b) who report chemical peculiarities in the spectrum 
of the HAeBe star W601 in NGC 6611. 

The only other star with a comparably young fractional 
age in which strong chemical peculiarities have been mea- 
sured in detail is the main sequence Bp star NGC 2244 334 
jBagnulo et all 120041 ). This star has been on the main se- 
quence for 2.3±0.3 Myr, but is more massive and hence prob- 
ably more evol ved than HD 72106A, with a fractional age 



= 0.02±0.01 jHensberge et al.ll2000l ; fBagnulo et ai1l2004 
Bagn ulo et al.1 (2004) find that the star has a temperature of 
15000 ± 1000 K and an observed longitudinal magnetic field 
strength of 9000 G. The authors find a strong underabun- 
dance of He by ~ 1 dex, as well as strong overabundances 
of Si and Fe by ~ 1 dex, and Ti and Cr by ~ 2 dex. These 
peculiarities are similar to those of HD 72106A, suggesting 
that there are important similarities between these stars. 

We find clear evidence for surface abundance inhomo- 
geneities in HD 72106A. Magnetic Doppler Imaging suggests 
that there is a large spot of overabundance in Si, Ti, Cr, and 
Fe near the positive magnetic pole. These inhomogeneities 
are similar to those seen in Ap/Bp stars, further support- 
ing a link between magnetic HAeBe stars and Ap/Bp stars. 
The Doppler reconstructions presented here are the earliest 
stage of intermediate-mass stellar evolution ever mapped. 

HD 72106B, in contrast to HD 72106A, appears to be 
chemically normal. While this is not the first binary system 
containing a ma gnetic chemically pe culiar and a chemically 
normal star (e.g. Carrier et al.ll2002j) . it does raise the ques- 
tion of how such systems are produced. Being a young bi- 
nary, it is likely that the stars of HD 72106 formed at the 
same time from approximately the same material. The tem- 
peratures of these stars differ by only 2250 ± 1120 K, and the 
vsini values of the stars differ by only ~ 10 kms -1 . How- 
ever, one star is magnetic and chemically peculiar while the 
other is not. This suggests that, whatever mechanism gave 
rise to the difference between the stars, it must be rather 
sensitive to the particulars of the star's initial conditions. 

It is instructive to compare the characteristics of HD 



72106A to those of the magnetic HAeBe stars HD 104237 
and HD 1900 73. HD 10423 7 was observed to possess a mag- 
netic field bvlDonati et al.l (Il997j ) and confirmed by iDonatH 
l|200d) and lAlecian et alj (|2008cl ). with a longitudinal field 
strength of ~ 50 G (Donati, private communication). HD 
1042 37 has a mass of about 2.3 Mp and an age of about 2 
Myr (|van den Ancker et al.lll998l ). lAcke fe Waelkensl (|2004l ) 
performed an abundance analysis of HD 104237 using equiv- 
alent widths, and found approximately solar abundances 
for a range of elements, including Si, Cr, and Fe. The 
sta r HD 190073 was r eported to possess a magnetic field 
by ICatala et al.1 l|2007i), with a longitudin al magnetic field 
strength of 74 ± 10 G. ICatala et~aiT(|2007l ) derive a mass of 
2.85 ±0.25 Mq and an age o f 1.2 ± 0.6 Myr (me a sured from 
the birth line) for this star. lAcke fe Waelkensl (2004) also 
studied the surface chemistry of this star and found roughly 
solar abundances. 

Thus it appears that the majority of known magnetic 
HAeBe stars are chemically normal, though some peculiar 
stars s eem to exist (such as NGC 6611 W601. lAlecian et al.l 
2008b). This is in contrast to main sequence Ap/Bp stars, in 
which magnetic fields are nearly always found with chemical 
peculiarities. More analysis of chemical abundances in mag- 
netic HAeBe stars must be performed, with an eye to iden- 
tifying chemically peculiar objects. Interestingly, the young 
chemically peculiar stars HD 72106A and NGC 2244 334 
display no emission, while HD 104237 and HD 190073 both 
display significant emission in their spectra. This suggests 
that HD 104237 and HD 190073 may still be undergoing sig- 
nificant accretion or mass loss while HD 72106A and NGC 
2244 334 are not. Thus it may be that accretion or mass 
loss mixes the stellar atmosphere, inhibiting the buildup of 
chemical peculiarities through diffusion. However, once ac- 
cretion halts chemical peculiarities may arise quickly. 



ACKNOWLEDGMENTS 

Thanks to Brian Mason at the United States Naval Obser- 
vatory for providing the Washington Double Star Catalogue 
data on HD 72106, as well as commentary on that data. 
CPF and CAW acknowledge support from the Academic Re- 
search Programme of the Canadian Department of National 
Defence. CPF, GAW, DH and JDL acknowledge Discovery 
Grant support from the Natural Sciences and Engineering 
Research Council of Canada. EA is supported by the Marie 
Curie FP6 program. 



REFERENCES 

Acke B., Waelkens C, 2004, A&A, 427, 1009 
Alecian E., Catala C, Wade G. A., Donati J.-F., Petit P., 
Landstreet J. D., Bohm T., Bouret J.-C, Bagnulo S., Fol- 
som C, Grunhut J., Silvester J., 2008a, MNRAS, 385, 
391 

Alecian E., Wade G. A., Catala C, Bagnulo S., Boehm T., 
Bohlender D., Bouret J.-C, Donati J.-F., Folsom C. P., 
Grunhut J., Landstreet J. D., 2008b, A&A, 481, L99 

Alecian E., Wade G. A., Catala C, Folsom C, Grunhut J., 
Donati J.-F., Petit P., Bagnulo S., Marsden S. C, Ramirez 
J., Landstreet J. D., Boehm T., Bouret J.-C, Silvester J., 



Magnetism and chemical peculiarities in HD 72106 15 



2008c, in Proceedings of the CP#AP Workshop Vol. 38 of 
Contributions of the Astronomical Observatory Skalnate 
Pleso, Magnetism in pre-MS intermediate-mass stars and 
the fossil field hypothesis, pp 235-244 
Alecian G., Stift M. J., 2007, A&A, 475, 659 
Bagnulo S., Hensberge H., Landstreet J. D., Szeifert T., 

Wade G. A., 2004, A&A, 416, 1149 
Carrier F., North P., Udry S., Babel J., 2002, A&A, 394, 
151 

Catala C, Alecian E., Donati J.-F., Wade G. A., Land- 
street J. D., Bohm T., Bourct J.-C, Bagnulo S., Folsom 
C, Silvester J., 2007, A&A, 462, 293 

Donati J.-F., 2000, These d'habilitation, Observatoire 
Midi-Pyrenees 

Donati J.-F., Mengel M., Carter B. D., Marsden S., Collier 

Cameron A., Wichmann R., 2000, MNRAS, 316, 699 
Donati J.-F., Semel M., Carter B. D., Rees D. E., Collier 

Cameron A., 1997, MNRAS, 291, 658 
Donati J.-F., Semel M., Rees D. E., 1992, A&A, 265, 669 
Donati J.-F., Wade G. A., Babel J., Henrichs H. f., de Jong 

J. A., Harries T. J., 2001, MNRAS, 326, 1265 
ESA, 1997, The Hipparcos and Tycho Catalogues, ESA 

SP-1200 

Fabricius C, Makarov V. V., 2000, A&A, 356, 141 

Gray D. F., 2005, The Observation and Analysis of Stel- 
lar Photospheres, 3rd edn. Cambridge University Press, 
Cambridge, UK 

Grevesse N., Asplund M., Sauval A. J., 2005, in Alecian 
G., Richard O., Vauclair S., eds, Element Stratification in 
Stars: 40 Years of Atomic Diffusion Vol. 17 of EAS Pub. 
Ser., The New Solar Chemical Composition, p. 21 

Hensberge H., Pavlovski K., Verschueren W., 2000, A&A, 
358, 553 

Hubrig S., Scholler M., Yudin R. V., 2004, A&A, 428, LI 

Jaschek C, Jaschek M., 1995, The behavior of chemical 
elements in stars. Cambridge University Press, Cambridge 

Kochukhov O., Piskunov N., 2002, A&A, 388, 868 

Kochukhov O., Piskunov N., Ilyin I., Ilyina S., Tuominen 
L, 2002, A&A, 389, 420 

Kupka F., Piskunov N., Ryabchikova T. A., Stempels H. C, 
Weiss W. W., 1999, A&AS, 138, 119 

Kurucz R., 1993, CDROM Model Distribution, Smithso- 
nian Astrophys. Obs. 

Landstreet J. D., Bagnulo S., Andretta V., Fossati L., Ma- 
son E., Silaj J., Wade G. A., 2007, A&A, 470, 685 

Lueftinger T., Kuschnig R., Piskunov N. E., Weiss W. W., 
2003, A&A, 406, 1033 

Michaud G., 1970, ApJ, 160, 641 

Michaud G., Char land Y., Megessier C, 1981, A&A, 103, 
244 

Morel P., 1997, A&AS, 124, 597 

Oudmaijer R. D., van der Veen W. E. C. J., Waters 

L. B. F. M., Trams N. R., Waelkens C, Engelsman E., 

1992, A&AS, 96, 625 
Palla F., Stahler S. W., 1993, ApJ, 418, 414 
Piskunov N., Kochukhov O., 2002, A&A, 381, 736 
Press W. H., Teukolsky S. A., Vetterling W. T., Flannery 

B. P., 1992, Numerical Recipes in FORTRAN, 2nd edn. 

Cambridge University Press, Cambridge 
Schiitz O., Meeus G., Sterzik M. F., 2005, A&A, 431, 165 
Shorlin S. L. S., Wade G. A., Donati J.-F., Landstreet J. D., 

Petit P., Sigut T. A. A., Strasser S., 2002, A&A, 392, 637 



Stibbs D. W. N., 1950, MNRAS, 110, 395 

Tikhonov A. N., 1963, Soviet Math. Dokl., 4, 1624 

Torres C. A. O., Quast G., de La Reza R., Gregorio-Hetem 
J., Lepine J. R. D., 1995, AJ, 109, 2146 

van den Ancker M. E., de Winter D., Tjin A Djie H. R. E., 
1998, A&A, 330, 145 

van Leeuwen F., 2007a, Hipparcos, the New Reduction of 
the Raw Data. Springer, Dordrecht 

van Leeuwen F., 2007b, A&A, 474, 653 

Vieira S. L. A., Corradi W. J. B., Alencar S. H. P., Mendes 
L. T. S., Torres C. A. O., Quast G. R., Guimaraes M. M., 
da Silva L., 2003, AJ, 126, 2971 

Wade G. A., Bagnulo S., Drouin D., Landstreet J. D., 
Monin D., 2007, MNRAS, 376, 1145 

Wade G. A., Donati J.-F., Landstreet J. D., Shorlin S. L. S., 
2000, MNRAS, 313, 851 

Wade G. A., Drouin D., Bagnulo S., Landstreet J. D., Ma- 
son E., Silvester J., Alecian E., B6hm T., Bouret J.-C, 
Catala C, Donati J.-F., 2005, A&A, 442, L31 



